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Abstract 

In this work, we investigate the possibility of enhancing the thermoelectric power (Seebeck coefficient) in graphene devices by 
strain and doping engineering. While a local strain can result in the misalignment of Dirac cones of different graphene sections 
in the ^-space, doping engineering leads to their displacement in energy. By combining these two effects, we demonstrate that a 
conduction gap as large as a few hundreds meV can be achieved and hence the enhanced Seebeck coefficient can reach a value 
higher than 1.4 mV/K in graphene doped heterojunctions with a locally strained area. Such hetero-channels appear to be very 
promising for enlarging the applications of graphene devices as in strain and thermal sensors. 
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1. Introduction 

The thermoelectric effect can be used to directly convert a 
temperature difference to an electric voltage and vice versa. 
When a conductor is connected to a hot and a cold reservoir 
with a temperature difference AT, an electrical voltage AV is 
established across the conductor according to 

AV = SAT (1) 


However, since the pioneering works of Hicks and Dressel- 
haus El, nanostructuring materials into low-dimensional sys¬ 
tems are now widely investigated to enhance the thermoelec¬ 
tric properties. To basically understand this size effect on the 
Seebeck coefficient, it is convenient to start from the simplified 
form of (2) derived for degenerately doped materials, i.e. 
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where S is the Seebeck coefficient characterizing the thermo¬ 
electric sensitivity of the conductor. The use of materials with 
high Seebeck coefficient is thus one of important factors to de¬ 
sign efficient thermoelectric generators and coolers or thermal 
sensors. It is also important to maximize the power factor crS^ 
where cr is the conductivity of the material. In electronic mate¬ 
rials in weak scattering regime, the linear response thermoelec¬ 
tric coefficient is given by the Mott’s formula d 
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where cr(6) is the energy-dependent conductivity associated to 
the density n(6) of electrons that fill energy states between e and 
evde, and /(e) is the Fermi-Dirac distribution function with the 
Fermi energy Ep. In conventional materials, a high Seebeck co¬ 
efficient is usually found in low carrier density semiconductors 
while a high conductivity is found in metals. The best com¬ 
promise is often to use heavily-doped semiconductors where, 
thanks to the finite bandgap, electrons and holes can be sepa¬ 
rated and the Seebeck coefficient is not reduced by their oppo¬ 
site contributions. 
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This expression suggests that any effect that can enhance the 
energy-dependence of the conductivity should enhance the See¬ 
beck coefficient, e.g., by enhancing the energy-dependence 
of the density n(6) that is directly dependent on the den¬ 
sity of states g(6). Hence, compared to bulk materials, low¬ 
dimensional systems are expected to provide higher Seebeck 
coefficient and power factor thanks to much higher dg(e)lde. 
For instance, it has been confirmed experimentally first in 
PbTe/Pbi-jcEujcTe quantum well structures m 

The potential of graphene as thermoelectric material is quite 
intriguing n This single layer of carbon atoms arranged in 
a honeycomb lattice offers fascinating electronic properties re¬ 
sulting in high mobility for massless chiral particles Ena. 
Regarding thermoelectric properties, graphene has the advan¬ 
tage of a strong energy-dependence of the conductivity near the 
charge neutrality point El However, it has the strong drawback 
to be gapless, which makes it difficult to separate the opposite 
contributions of electrons and holes to the Seebeck coefficient. 
It results in a finite but small value of A < 100 juV/K in pristine 
graphene lO. 

So far, many studies have suggested different ways to open 
a band gap in graphene. As a direct consequence, it has been 
shown that the Seebeck effect can be significantly enhanced in 
graphene nanostructures with finite energy gap such as arm¬ 
chair graphene nanoribbons IQ, hybrid structures combining 
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zigzag graphene nanoribbon with zigzag boron nitride nanorib¬ 
bon IITQI . graphene nano-hole lattices ITll . graphene nanorib¬ 
bons consisting of alternate zigzag and armchair sections ifT^ . 
vertical graphene junctions C3, and graphene p-n junctions 
na. However, each mehod has its own drawbacks and still 
need to be confirmed by experiments. 

Furthermore, graphene has been demonstrated to be con¬ 
formable and able to sustain large strain GSHnl, making it a 
promising candidate for flexible devices. Also, strain engineer¬ 
ing has been proposed to be an alternative approach to modulat¬ 
ing the electronic properties of this material. In particular, it has 
been shown that a gap can be opened in pristine graphene for 
deformations beyond 20 % (m. In recent works (1911201, we 
have investigated the effects of uniaxial strain on the transport 
properties of 2D graphene heterochannels and found that a sig¬ 
nificant conduction gap of a few hundred meV can be achieved 
with a small strain of a few percent. This conduction gap is 
not due to a bandgap opening in the band structure but to the 
strain-induced shift of the Dirac cones in the Brillouin zone of 
different graphene sections. This result motivated us to investi¬ 
gate here the possible strain-induced enhancement of Seebeck 
coefficient in graphene nanostructures. In addition, doping en¬ 
gineering has been included in our investigation since it is likely 
to increase strongly the conduction gap, and thus the Seebeck 
coefficient, in graphene doped heterojunctions. 

Regarding some thermoelectric applications, the figure of 
merit ZT is another important parameter. It is defined as 
Zr = crS^T Ik, where k is the thermal conductivity. Actually, 
the strain engineering is not an effecient technique to modulate 
the phonon bandstructure ED and to strongly reduce the ther¬ 
mal conductivity in the junctions studied here. Hence, though 
the Seebeck coefficient and the power factor crS ^ are strongly 



Figure 1: Schematic view of the graphene structure investigated in this work, 
consisting of a PN diode with a locally strained area of length Ls that covers 
symmetrically both doped sides. The bottom shows its strain and potential 
profiles where the doping is characterized by the potential difference At/ = 
Ur - Ul and the length Lj of transition region. 


improved, we believe that the combination of this design with 
additional nanostructuring (e.g., as in ifTSll ) or more complex 
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Figure 2: (a) conductance and (b) Seebeck coefficient as a function of Fermi 
energy Ep for different strain amplitudes in the device with uniform doping, 
(c) maximum value of Seebeck coefficient (S max) and conductance (Goff) in 
OFF state (i.e. at = 0) as a function of strain amplitude. 
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design would be required to achieve high ZT. For this reason, 
we focus here our investigation on the Seebeck coefficient that 
is an essential ingredient. 


2. Model and calculations 

We investigate 2D graphene doped heterojunctions with a 
strain area of finite length as schematized in Fig. 1. The strain 
area covers symmetrically both doped sides and its length Lg 
is assumed to be much longer than the length of the transi¬ 
tion region Lj between left and right doped sections. It has 
been shown that the doping profile can be generated/controlled 
by chemical doping or electrostatic methods, e.g., see refs. 
(^1231. Though expected to be short for achieving high band- 
to-band tunneling current, the transition length is always finite 
in the devices with chemical doping Ea. In the case of electro¬ 
static doping (231, this length is also finite but can be controlled 
by tuning the properties of insulator layer, i.e., its thickness 
and dielectric constant. Throughout this work, unless otherwise 
stated our calculations were performed at room temperature for 
Lg = 10 nm and Lj = 10 nm. 

A -orbital tight-binding model was used to calculate elec¬ 
tronic and thermoelectric properties of the device. The Hamil¬ 
tonian is Htb = l^nmtnmclcm where tnm is the hopping energy 
between nearest neighbor atoms. We consider a local uniaxial 
strain applied along the Oy direction. Accordingly, the strain- 
dependence of C - C bond vectors is given by 

/ = (1 - w)^x(O) ... 

\ r,(^) = (l+^)r,(0) 

where r] is the strain amplitude and y = 0.165 is the Poisson’s 
ratio (241 . The hopping interaction between atoms is defined by 
tnmili) = toexp[-337(rnm(rj)lro - 1)] (HI, where = -2.1 eV 
and TnmiO) = vq = 0.142 nm are the hopping energy and the C-C 
distance in the unstrained case, respectively. In the specific case 
of this study, where the strain is applied along a zigzag axis, we 
have to consider two different hopping parameters 0,2 in the 
armchair and zigzag directions, respectively, corresponding to 
two bond vectors ri^ in strained graphene as in (251 . 

This tight-binding Hamiltonian was solved by 
the Green’s function technique (^ : Q{e,ky^ = 

[e + - Htb (e, ky) - Er (e, ky)j where H,b (ky) is 

the Hamiltonian rewritten in the wavevector /: 3 ;-dependent 
(quasi-ID) form as in (T9l and ^l(R) is the self en¬ 
ergy describing the left (right) contact-to-device cou¬ 
pling. The local density of states and the transmission 
probability needed to evaluate the transport quantities 
are determined as D(^€,ky,rn) = jn 

and Te{e,ky^ = Tr |Fl^F/?^’^|, respectively, where 

^L{R) = i{^L(R) is the transfer rate at the left (right) 

contact. The electrical conductance and the Seebeck coefficient 
were calculated by 


G(Ef) = GoLo(Ef,T) (5) 



- 0.6 - 0.4 - 0.2 0 0.2 0.4 0.6 


Fermi energy (eV) 



Fermi energy (eV) 



Figure 3: (a) conductance and (b) Seebeck coefficient as a function of Fermi 
energy Ep for different AU. (c) maximum Seebeck coefficient (Smax) and con¬ 
duction gap (Egap) as a function of AU. Inset: Smax as a function of Egap. 
77 = 5 % is considered here. 
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Figure 4: Local density of states in the devices of different At/: 0 (a), 0.2 eY (b), and 0.6 eV (c). The blue color regions correspond to energy-gaps, i.e., low density 
of states, (d) shows the transmission coefficient in these three cases, ky - and 77 = 5 % are considered here. 


where 

, T)=^J dKydeTeie, Ky){e - Epf j. (7) 

Here, Go = e^WjhLy and the channel width W = MceiiLy with 
the size of unit cells Ly (= tq ^^3) and the number of cells Mceii 
along Oy direction. The integral over Ky (Ky = kyLy) is per¬ 
formed in the whole Brillouin zone. The distribution function 
/ ( 6 , Ep) = (1 + exp [(6 - Ep)/kBT])~^ is the Fermi-Dirac func¬ 
tion with the Fermi energy Ep. 

3. Results and discussion 

3.1. Strained device with uniform doping 

First, let us examine the basic effect of strain on the transport 
properties of graphene. Actually, for a small strain of a few per¬ 
cent, graphene is still metallic IT^ . i.e., its gapless character do 
not change. However, even a small strain causes a shift of Dirac 
points in the space ca. As a consequence, it may lead to the 
opening of a conduction gap in strained/unstrained graphene 
junctions mu. This phenomenon is explained as follows. For 
a given ^^-mode, strained graphene has a ^^-dependent energy 

gap expressed as E^^^^iky) = 2 | 2/2 |cos(^)| - /i|. In the un¬ 
strained graphene, E^^^^^^^iky) has a similar expression but with 

= /2 = ^ 0 - Due to these effects of strain, the Dirac cones 
of unstrained/strained graphene are formed at different posi¬ 
tions in the ^^-axis and hence a finite energy gap E^^^iky) of 
the device transmission is always obtained for all wavevectors 
ky. The value of E^^Jp(ky) corresponds to the maximum value of 
^u^stramQ^^'^ and Eg^p^(ky). Finally, the conduction gap is given 
by the minimum value of Eg^^p(ky) in the whole Brillouin zone, 
i.e, Egap = 2|^^/o| in the present case 1^ . Besides, it has 


been also shown that the properties of conduction gap strongly 
depend on the amplitude of strain, its applied direction and lat¬ 
tice orientation csi. 

In Figs. 2(a) and 2(b), we plot the conductance and Seebeck 
coefficient, respectively, as a function of Fermi energy for dif¬ 
ferent strain amplitudes p ranging from 0 to 6 %. While the 
minimum value of conductance is finite in pristine graphene 
(77 = 0 ), a finite conduction gap is achieved when a local strain 
is applied to the structure. Actually, Egap increases almost lin¬ 
early as a function of strain amplitude and, particularly, it takes 
the value of 0.162, 0.324 and 0.486 eV for 77 = 2 %, 4 %, and 
6 %, respectively. As expected from this gap enhancement, 
the maximum value S max of Seebeck coefficient increases from 
0.086 mV/K for 77 = 0, in agreement with experimental data El, 
to 0.803 mV/K for 77 = 6 %. We find as shown in Fig. 2(c) that 
the conductance in the OFF state, i.e., the minimum of conduc¬ 
tance when varying Ep (practically, dXEp = 0), decreases ex¬ 
ponentially when the strain amplitude increases, in accordance 
with the linear increase of conduction gap. This effect is at the 
origin of the enhancement of ON/OFF current ratio in the tran¬ 
sistors based on this type of strain heterochannels (251, where 
the OFF and ON currents are the minimum and maximum cur¬ 
rents, respectively, obtained when tuning the gate voltage. It is 
shown concomitantly in Fig. 2 (c) that S max increases linearly 
as a function of strain amplitude and reaches up to 1.35 mV/K 
for 77 = 10%. However, practically it may be difficult to ap¬ 
ply locally such a large strain in this type of structure. Hence, 
we propose here to introduce appropriate doping engineering 
in this device to enhance the conduction gap with a reasonable 
strain amplitude. 

3.2. Devices with both strain and doping engineering 

In this subsection, we discuss the effects of doping engi¬ 
neering schematized in Fig. 1 on the transport properties of 
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this strain heterostructure. The structure now consists of three 
sections: different doped-graphene sections in both sides and 
a transition region between them. Note that within the con¬ 
dition Lg » Lt, left and right doped sections are actually 
formed by two different strain junctions. In Fig. 3, we display 
the conductance G and the Seebeck coefficient as a function of 
Fermi energy Ep for different doping profiles. The strain am¬ 
plitude is fixed to 5 % and the doping profile is characterized 
by a finite potential difference AU = Ur - Up (see Fig. 1). 
The doping engineering consists here in controlling the carrier 
density profile, which is characterized and determined by both 
the potential AU and the Fermi level Ep. The results in Fig. 
3(a) show that for small AU, the conduction gap increases with 
this potential difference, which can be explained as follows. 
When AU increases, Egap of strained/unstrained junctions in 
the left and right sides are shifted in opposite directions, which 
results in the enlargement of the transmission gap (i.e., con- 




Figure 5: (a) conductance and (b) Seebeck coefficient as a function of Fermi 
energy Ep with high values of AC/. 77 = 5 % is considered here. 


duction gap). Indeed, this phenomenon is clearly illustrated in 
the pictures of local density of states (LDOS) and transmission 
coefficient in Fig. 4, obtained for ky = (K^^strain ^ K^train^^i2 

where strain wavevector at the Dirac points of 

unstrained/strained graphene, respectively. Note that at this ky 
point, E-^J^^^^ky) - - Egap for AU = 0. Actually, 

the conduction gap Egap{AU) of the whole device in the range 
of small AU (i.e., AU < EgapiU)) is basically determined as 
Egap{AU) = Egap{U) + AU, as shown in Fig. 4(d) and confirmed 
in Fig. 3(c) where we find the linear dependence of Egap as a 
function of At/. As a consequence, S is significantly enhanced 
when increasing AU (see Fig. 3(b)) and Smax reaches the value 
of 0.67, 0.84, 1.0, 1.16 and 1.32 mV/K for AU = 0, 0.1, 0.2, 
0.3 and 0.4 eV, respectively. It is worth noting here that the 
result obtained for rj = 5 %, AU = 0.4 is almost similar to 
that obtained forr]= 10%, AU = 0, which is about 15 times 
greater than Smax in pristine graphene. Thus, it is demonstrated 
that the doping engineering can be an effective way to further 
enlarge Egap without the requirement of too large strain. All the 
features above are clearly summarized in Fig. 3(c). In addition, 
the inset of Fig. 3(c) confirms that whatever the technique of 
energy gap opening Smax always depends linearly on Egap, as 
predicted theoretically in ref. ifTOl . 

Next, we go to analyze the effects of large AU, i.e., values 
greater than EgapiO). In Figs. 5(a) and 5(b) we plot the conduc¬ 
tance and Seebeck coefficient as a function oiEp for large val¬ 
ues of AU increasing from 0.4 eV to 0.6 eV. It is shown that the 
conduction gap is separated in two smaller ones that correspond 
to the conduction gap of each strained/unstrained junction of the 
structure. Between these two gaps, a region of finite conduc¬ 
tance is recovered due to the band-to-band tunneling (BTBT), 
as in a standard doped tunnel diode, which is illustrated clearly 
in Figs. 4(c) and 4(d) where we plot the map of LDOS and 
transmission coefficient, respectively, for AU = 0.6 ^V. As a 
consequence, at large AU the Seebeck coefficient exhibits two 
positive (negative) peaks with a maximum value Smax that re¬ 
duces and tends finally to a finite value Soo when increasing 
AU, as shown in the inset of Fig. 5(b). Note that the value Soo 
is generally higher than the value of Smax obtained for Af/ = 0 
but tends to this value when the transition length Lp is reduced. 
This is explained by the detrimental contribution of band-to- 
band tunneling current, which is significantly reduced when 
increasing Lp (see further comments below). It is addition¬ 
ally shown that Smaxi^U) has a peak at Af/ = EgapiO), e.g., 
AU ^ 0.4 eV for 7] = 5% here. 

Finally, we would like to clarify the roles of device parame¬ 
ters Eg and Lp on the obtained results. Note that in this work, 
we consider only the case Eg » Lp. Within this condition, 
each doped section contains by a strain junction and hence has 
a finite conduction gap. In principle, this conduction gap is 
strongly dependent on the length of the strained graphene part 
in these two doped sections, i.e., the transmission probability 
in the gap increases exponentially when reducing the length of 
strained graphene area. Hence, to ensure that the transmission 
is fully suppressed in the gap, the length Eg should be much 
larger than Lp. Moreover, the required length Eg is basically 
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dependent on the value of energy gap, i.e., the larger Ls is re¬ 
quired for the smaller Egap (i.e., smaller strain). In particular, 
in the case of 77 = 5 %, > Lj + 20 nm should be used. Ad¬ 

ditionally, it has been shown that the length Lp of the transi¬ 
tion region between n- and p-doped sections plays an impor¬ 
tant role on the BTBT current 1271 . i.e., this current is exponen¬ 
tially reduced when increasing Lp as seen in Fig. 6 (a). More 
interestingly, based on this reduction of BTBT current, the See- 
beck coefficient in the case of Af/ > EgapiO) is significantly 
enhanced when increasing Lp, i.e., Smax reaches 1.48 mV/K for 
Lp = 40 nm as shown in Fig. 6 (b) while it is only about 0.77 
mV/K for Lp = 5 nm. We notice that in the case of small AU, 
the BTBT current is negligible as shown for At/ = 0 and 0.2 
in Fig. 4 and hence S is very weakly dependent on Lp. 



Fermi energy E|_ (eV) 



Figure 6: (a) conductance and (b) Seebeck coefficient as a function of Fermi 
energy Ep for the different lengths Lp. rj = 5% and AU = 0.5 eV. 


4. Conclusion 

In this work, we have proposed to make appropriate use of 
strain and doping engineering to generate and enlarge a con¬ 
duction gap in graphene heterochannels and to benefit from this 
feature to enhance the Seebeck effect. The maximum value 
S max of the Seebeck coefficient was shown to increase linearly 
with the conduction gap. Remarkably, with a small strain of 
5 % and an appropriate doping profile, the Seebeck coefficient 
can reach a value higher than 1.4 mV/K, i.e., 17 times higher 
than the value in gapless pristine graphene. Besides its use in 
strain sensors, this design strategy is thus promising to achieve 
good performance in graphene devices based on the Seebeck 
effect, as thermal sensors. 
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